Degenerate three-level laser with parametric amplifier and squeezed vacuum 
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Applying stochastic differential equations, we study the squeezing and statistical properties of 
the cavity and output modes of a degenerate three-level laser whose cavity contains a parametric 
amplifier and coupled to a squeezed vacuum reservoir. We consider the case in which the top 
and bottom levels of the three-level cascade atoms injected into the cavity are coupled by the pump 
mode emerging from the parametric amplifier. It turns out that the presence of the squeezed vacuum 
reservoir and the parametric amplifier contribute considerably to the mean photon number and the 
degree of squeezing of the cavity and output modes. It appears that almost perfect squeezing can 
be achieved at steady state and at threshold for a suitable choice of parameters. 
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I. INTRODUCTION 



The quantum properties of three-level lasers have been investigated by several authors [l|, 0, 0, 0, 0, 0, 0, 0, 0, El • 
It is found that three- level lasers can generate squeezed light under certain conditions. The generation of squeezed 
^ ' light by three-level lasers could be realized when either the atoms are initially prepared in a coherent superposition 
\0 • of the top and bottom levels [l|, 0, 0, 0, 0, 0, EU or when these levels are coupled by a strong coherent light 0, 0] . 

More recently a three-level laser with a parametric amplifier has been studied when either the three-level atoms are 
initially prepared in a coherent superposition of the top and bottom levels [l[ or when these levels are coupled by the 
pump mode emerging from the parametric amplifier Q . These studies show that the effect of the parametric amplifier 
is to increase both the mean photon number and the intracavity as well as the output mode squeezing significantly. 
All previous studies have been confined to the case for which the cavity mode of the three-level laser is coupled to a 
vacuum reservoir. Moreover, apart from the squeezing spectrum of the output mode, all attention has been paid to 
the calculation of the squeezing and statistical properties of the cavity mode. Since the output mode is accessible to 



I ' measurement, it appears to be appropriate to study the quantum properties of this mode. 



In this paper we analyze the squeezing and statistical properties of the cavity and output modes of a degenerate 
three-level laser whose cavity contains a degenerate parametric amplifier and coupled to a squeezed vacuum reservoir 
via a single-port mirror. We consider the case for which the top and bottom levels of the three-level atoms injected into 
CJ^ 1 the cavity are coupled by the pump mode emerging from the parametric amplifier. We obtain stochastic differential 
| equations for the cavity mode variables employing the pertinent master equation. Using the solutions of these 
equations, we calculate the quadrature variance, the mean photon number, and the power spectrum of the cavity and 
output modes. In addition, we determine the squeezing spectrum of the output mode. 



II. STOCHASTIC DIFFERENTIAL EQUATIONS 



We consider a laser cavity containing a parametric amplifier and coupled to a squeezed vacuum reservoir. Three- 
level atoms in a cascade configuration are injected into the cavity at some constant rate r a . We denote the upper 
level by \a), the middle level by \b), and the lower level by |c), as shown in Fig. 1. The diploe allowed transitions 
between levels \a) and \b) and between levels \b) and |c) are resonant with the cavity mode. The direct transitions 
between levels \a) and |c) are dipole forbidden. 

With the pump mode treated classically, a degenerate parametric amplifier is describable in the interaction picture 
by the Hamiltonian 

H= l ±(aV-a% (1) 
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in which e = Xp with A and p being respectively the coupling constant and the amplitude of the pump mode. The 
equation of evolution of the density operator associated with this Hamiltonian has the form 



= |(/5a 2 - a 2 p + a) 2 p - pa) 2 ). 



(2) 



In addition, the Hamiltonian describing the interaction of a three-level atom with the cavity mode and with the pump 
mode emerging from the parametric amplifier has the form 



H = ig[tf(\b)(a\ + \c)(b\) a(\a){b\ + \b)(c\)] + &\c)(a\ - \a)(c\) 



(3) 



where is proportional to the amplitude of the pump mode and g is the atom-cavity mode coupling constant. We 
take the initial state of a single three-level atom to be 



|Vu(0)) = -L|a)+i=|c) 



and hence the density operator for a single atom is 



^(0) = i|a)<a| + i|a>(c| + i| C )< a | + i|c)<c|. 



(4) 



(5) 



Following the procedure developed in Ref. Q , we can show that the master equation for the laser cavity mode coupled 
to a squeezed vacuum reservoir to be 



dt 



p = A(2a) pa — aa) p — paa)) + B(2apa) — a) ap — pa) a) 
+ C(a)pa) + apa — pa) 2 — a 2 p) + T>(a) pa) 4- apa — pa 2 — a) 2 p), 



(6) 



where 
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kN A 
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2 + 4B 



3/3 



(7b) 



kM A 
~2~ + 4B 



2 2 
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(7c) 



V 



kM A 
~ + 4B 



(3 (3 2 (3 3 



(7d) 



B = (l+/3 2 )(1 + ^), 



N = sinh r, 



M = sinh r cosh r, 



(7e) 
(7f) 

(7g) 
(7h) 



(8) 
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is the linear gain coefficient, k is the cavity damping constant, 7 is the atomic decay rate assumed to be the same for 
all the three levels, and r is the squeeze parameter. Therefore, on account of Eqs. and © the master equation 
for the cavity mode of the quantum optical system under consideration becomes 

— p = ^P^ 2 ~ ^ 2 P ~ P^ 2 ) + A(2a t pa — da) p — pad)) + B{2dpd) — d) dp — pd) d) 

+ C(d) pd) + apa — pa) 2 — a 2 p) + T>(d) pa) + dpd — pa 2 — a) 2 p). (9) 

We next seek to determine, applying this master equation, stochastic differential equations for the cavity mode 
variables. To this end, applying ([9]) one readily finds 

±{a) = -{B-A){a) + {C-V + e){a)), (10) 
at 

4-(a 2 ) = -2{B - A)(a 2 ) + 2(C - V + e)(a)a) + e - 22?, (11) 
at 

j{a)a) = -2(B - A){tf&) + {C-V + e)«a t2 ) + (a 2 )) + 2A (12) 
and the c-number equations corresponding to these normally ordered equations are 

±{a) = -(B-A){a) + (C -V + e){a*), (13) 
at 

4(a 2 ) = -2(B - A)(a 2 ) + 2(C - V + e)(a*a) +e-2V, (14) 
at 

j t (a*a) = -2{B - A)(a*a) + {C-V + e){(a* 2 ) + (a 2 )) + 2A. (15) 
Based on Eq. (|13[) . one can write 



^-a(t) = -(B - A)a{t) + (C-V + e)a*{t) + f{t), (16) 
at 

where f(t) is a noise force the properties of which remain to be determined. We note that Eq. (fl3|) and the expectation 
value of Eq. (|16| will have the same form provided that 



(/(*)) = 0. (17) 

It is easy to see that 

±(a 2 (t))=2(a(t)±a(t)) (18) 



and on substituting Eq. (fl6|) into (jTSJ) , we get 

d 



^<a 2 (i)) = -2(6 - A)(a 2 (t)) + 2(a(t)f(t)) +2(C-V + s)(a*(t)a(t)). (19) 

It can also be verified in a similar manner that 

j t {a*(t)a(t)) = -2(B A){a*(t)a(t)) +(C-V + e)({a* 2 (t)) + {a 2 (t))) 

+ (a(t)r(t)) + (a*(t)f(t)). (20) 
Comparison of Eqs. (TT4")) and (fT5]) as well as Eqs. (TT5|) and (T2T))) leads 

(a(t)f(t)) = ±(e-ZD), (21) 



(a(t)f* (t)> + («*/(*)) = 2A (22) 
Furthermore, one can write a formal solution of Eq. (|16[) as 



a(t) = a{Q)e- {B - A)t + / e -(S-^)(t-t') [( C _ p + e ) a *( t ') + f(t')]dt'. (23) 
j o 

Multiplying Eq. (f2"3"|) on the left by /(i) and taking the expectation value, we have 

(«(*)/(*)) - (a(0)/(i))e-( B - A )' + /' e -( B -^(*-*')[(C - P + e)<a*(i')/(i)) + (f{t')f{t))]dt'. (24) 



t 

Assuming that the noise force at time t does not affect the cavity mode variables at earlier times and on taking into 
account (|2T|) . Eq. (|24| can be put in the form 

* e -(B-A)(t-t>) {f(j ., )f ( t))dtl = 1 (£ _ ^ (25) 
^ 

On the basis of this result, one can write @ 

(f(t')f(t)) = (e-2V)5(t-t'). (26) 

It can also be established in a similar manner that 

(f(t)r(t>))=2A6(t-t>). (27) 

We note that Eqs. (|26|) and ([27)1 describe the correlation properties of the noise force f(t) associated with the normal 
ordering. 

Now introducing a new variable defined by 

a±(t) = a*(t)±a(t), (28) 

we easily get with the help of (16) that 

d_ 

dt 

where 



a±(t) = -A T « ± (() + f(f)±/(t), (29) 



\ T = (B-A)T(C-T> + e). (30) 

The solution of Eq. (f2"9")l can be written as 



a±( 



(t) = a ± (0)e- x ^ + / e-^-t'Hrit') ± f(t'))dt'. (31) 



It then follows that 



in which 



and 



with 



a(t) = E+(t)a{0) + E_{t)a*(0) + F(t), (32a) 



E ± {t)= l -{e- x - t ±e- x + t ), (32b) 



F(t) =F+(t)+F_(t), (33a) 



F±(i) = 1 -J\- X ^ t ~ t '\r{t')±f(t'))dt'. (33b) 
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III. QUADRATURE VARIANCE 

In this section we seek to analyze the quadrature variance of the cavity and output modes. 

A. Quadrature variance of the cavity mode 

We define the quadrature operators for the cavity mode as 

a + = a) + a (34) 

and 

a_ = i{a) — a). (35) 

These quadrature operators satisfy the commutation relation [a+, a_] = 2i. The variance of these quadrature operators 
is expressible in terms of c-number variables associated with the normal ordering as 

Aa 2 ± = l±(a±(t),a±(t)), (36) 

where a± (t) is given by Eq. (|28[) ■ Assuming the cavity mode to be initially in a vacuum state and taking into account 
(|3"Tj) along with (fl~7|) . wee see that 



(«±(t))=0. (37) 

Thus in view of this result, Eq. ([36| reduces to 

Aa 2 ± = 1 ± (a<±(t)). (38) 

Furthermore, applying Eq. ([29]), one easily gets 

|<<4(t)> = -2A T <<4 (*)> + 2<a±(t)/*(t)) ± 2<a±(t)/(t)>. (39) 



With the aid of Eq. (J3Tj) along with (J26J) and (|27j). we readily obtain 



^ (a 2 ± (t)) = -2A T <a 2 ± (t)) + 2( £ - 223 ± 2^1) (40) 



and at steady state we have 



<<4(t)>„ = -. (41) 

Hence on account of Eqs. (|4"Tj) , (|5D|). and ([7a|) - (|7h|) . the quadrature variance takes at steady state the form 

2 = 2 K (l + /3 2 )(l + /3 2 /4)e 2 - + A(4 + /3 2 ) 

+ 2( K -2e)(l + /5 2 )(l + /5 2 /4)+A(2/3-/? 3 ) 1 j 

and 

2 = 2 K (l+/? 2 )(l + /? 2 /4) e - 2 '- + 3A/? 2 

- 2( K + 2 £ )(l + /3 2 )(l + /3 2 /4) + A(4 / 9 + /33)- 1 j 

We next proceed to calculate the quadrature variance of the cavity mode at threshold. We note that Eq. (f2"9"| will 
not have a well-behaved solution if A_ < 0. Hence we identify A- = as the threshold condition. With the aid of 
Eq. (f3"0"|) together with (fTa"|) - ([7h")) , the threshold condition can be written as 



l ^ h * _ (A1\ 

e = o + —i , Wlx ,2,. ' ( 43 ) 



k A(2/3 ~ P 3 

2 + 4(l + /3 2 )(l + /3 2 /4) ■ 

Now introducing this into Eqs. (I42ap and (|42b>[) . the quadrature variance at threshold takes the form 

Aa 2 + -> oo (44a) 
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and 



and 

a 2 = 2 K (l + /3 2 )(l + /3 2 /4) e - 2 '- + 3A/3 2 

4 K (l + /3 2 )(l+/? 2 /4) + 6A/3 ' 1 1 

Fig. 2 indicates that the quadrature variance increases with j3 and decreases with the squeeze parameter r. For 
r = 1, K = 0.8, and A — 100, the minimum value of the quadrature variance is found to be 0.022 at (3 — 0.022. We 
immediately see that the intracavity squeezing is 97.8% below the vacuum level. We also note from Fig. 3 that the 
presence of both the parametric amplifier and the squeezed vacuum reservoir increases the intracavity squeezing over 
and above the squeezing achievable due to the coherently driven three-level laser Q. Moreover, Fig. 3 clearly shows 
that the effect of the squeezed vacuum reservoir is significant for relatively small values of p. 

In the absence of the nonlinear crystal (NLC) the quantum optical system under consideration reduces to a 
coherently driven degenerate three-level laser coupled to a squeezed vacuum reservoir. The quadrature variance for 
this system takes upon setting e = A/i = (with [i ^ 0) in (|42a|) and (|42b[) the form 

2 2 K (l + /? 2 )(l+/? 2 /4)e 2 '- + A(4 + /3 2 ) 
°+ 2k(1 + /? 2 )(1+/? 2 /4) + A(2/3-/33) { ' 

2 2 K (l+/? 2 )(l+/? 2 /4) e - 2 '- + 3A/? 2 

- 2 K (l + /3 2 )(l + /3 2 /4)+A(4/3 + /33)- 1 1 

As can be seen from Fig. 4 the degree of squeezing increases with the squeeze parameter r and decreases with (3. Fig. 
5 also indicates that the presence of the squeezed vacuum reservoir increases the intracavity squeezing significantly 
for relatively small values of (3. For A=100, k = 0.8, and r = 1.0, the minimum value of the quadrature variance given 
by (|45b[) turns out to be 0.035 and occurs at p = 0.023. We then note from this result that the intracavity squeezing 
is 96.5%. 

B. Quadrature variance of the output mode 

The variance of the quadrature operators for the output mode, defined by 

a+.out = a\ ut + a out (46) 

and 

a-,out = i{a\ ut - a out ), (47) 
can be expressed in terms of c-number variables associated with the normal ordering as 

Aa± out = 1 ± (a± t0ut (t), a±,out(*)>i (48) 

in which 

a±,out(t) = y/K,a±(t) - a±,i n (t), (49) 

and 

a± liB (t) = 4=(/fl(*)± /«(*)), (50) 

with fn being the noise force associated with the squeezed vacuum reservoir. This noise force has the following 
correlation properties: 

(fn(t)) = 0, (51a) 

(mtW)} = KNS(t-t'), (51b) 

(f* R (t)m')) - (Mt)f R (t')) = K MS(t - t'). (51c) 
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In view of Eqs. (150)) and (I51a[) . one can write 

(a±, in (t)) = 0. 

Thus on account of Eqs. (|37|) and (l52|) . the quadrature variance (|48|) takes the form 



l ±,out 



= l±(a 



(*)>• 



Furthermore, with the aid of Eq. 



we find 



(a 2 ±tOUt (t)} = n{a 2 ± {t)) - 2^{a±(t)a±, in (t)) + (al in (t)). 



It can be readily verified that 



and 



(a ± (t)a±, in (t)) = \/k(M ± N) 



(al in (t))=2(M±N). 



Upon substituting Eqs. (gT]), (J55J), and ((56j) into (54)), we get 

(«±,„„tW) S s = K 



Aq 



+ 2(1-k)(M±7V). 



(52) 
(53) 
(54) 
(55) 
(56) 

(57) 



Therefore, with the aid of this result and along with (fTa|) - (|7h[) , the quadrature variance takes at steady state the form 
"2(1 + (3 2 ){l + 2 /A)[2e + n{e 2r - 1)] + A(4 - 2/3 + 1 + /3 3 )" 



,out 



= 1 + K 



and 



Aa-.out = 1 - « 



2(k - 2e)(l + /? 2 )(1 + /? 2 /4) + A(2/? - /3 3 ) 
2(1 + /3 2 )(1 + /3 2 /4)[2e + K (l - e - 2r )] + A(40 - 3/3 2 



/3 3 ) 



2(k + 2e)(1 + /3 2 )(1 + /3 2 /4) + A(4/3 + /3 3 ) 
On account of Eq. (|4"3"| the quadrature variance reduces at threshold to 



+ (l-«)(e*-l) 



(1-«)(1 



and 



1-K 



2(1 + /3 2 )(1 + /3 2 /4)(2 - er 2r ) + A{60 - 3/3 2 ) 



4k(1 + /3 2 )(1 + /3 2 /4) + 6A/3 



(l-«)(l-e- 2r ). 



(58a) 



(58b) 



(59a) 



(59b) 



We note from Fig. 6 that the amount of squeezing of the output mode is less than that of the cavity mode for 
relatively small values of (3. The minimum value of the quadrature variance of the output mode is 0.045 at j3 = 0.022. 
This shows that the maximum squeezing of the output mode is 95.5%. 

We next wish to consider the case for which the nonlinear crystal is removed from the cavity. To this end, upon 
setting e — A/i = (with fi ^ 0) in Eqs. (|58all and (|58bj) . the quadrature variance for this case reduces to 



A a+, «t = 1 + « 



2k(1 + /3 2 )(1 + /3 2 /A)(e 2r - 1) + A(4 -2f3 + f3 2 + (3 3 ) 



and 



Aa 2 



,out 



2k(1 + /3 2 )(1 + /3 2 /4) + A(2p - (3 3 ) 
2k(1 + (3 2 ){1 + /3 2 /4)(l - e- 2r ) + A(A/3 - 3/3 2 + f3 3 ) 



+ (l-n)(e 2r -l) 



2k(1 + (3 2 ){l + /3 2 /4) + A(4/3 + /3 3 ) 



(1-«)(1 



-2r\ 



(60a) 



(60b) 



We notice that Eqs. (|60a[) and (|60b[) represent the quadrature variance of the output mode of a coherently driven 
threedevel laser coupled to a squeezed vacuum reservoir. In Fig. 7, we plot Eqs. ((45b)) and ((60b)) versus /3. These 
plots also show that the degree of squeezing of the output mode is less than that of the cavity mode. We have found 
that the minimum value of the quadrature variance of the output mode to be 0.055 at j3 = 0.023. This indicates that 
the maximum squeezing of the output mode for this case is 94.5% below the vacuum level. We observe that the degree 
of squeezing of the of the output mode in the presence of the NLC is greater by 1% than that without the NLC. 



IV. SQUEEZING SPECTRUM 



In this section we calculate the squeezing spectrum of the output mode. The squeezing spectrum of a single-mode 
light is expressible in terms of c-number variables associated with the normal ordering as 



/>oo 

S£*(w) = l±2Re / « ut (i),<4 ut (t + T)) ss e lu)T dr. 
Jo 

Thus on account of Eqs. (|37p and (|52[) . the squeezing spectrum can be put in the form 

/>oo 

S^iuj) = l±2Re / (c4 ut (*)c4"*(f + T)) ss e iuJT dT, 
Jo 

so that introducing Eq. (|49|) . we get 

POC POO 

S^ ut (to) = 1 ± 2 K Re / (a±(t)a±(t + T)) ss e iu}T dr =p 2^i?e / (a±(i)c4*(* + r)) ss e MT dr 



'0 

fOC 



T2v^i?e / (a! t ™(i)Q !± (t + r)) ;s;j e^ r dT±2i?e / (a£(i)a± (i + r)) s;j e^ r dr. 
Jo Jo 



Furthermore, the solution of Eq. (|29|) can also be written as 



«±(i + t) = a±{t)e~ x ^ T + [ , 

Jo 



-X t (t-t') ( f* 



(r(t + T')±f(t + T'))dr'. 



Now applying this equation along with the quantum regression theorem, one can establish that 

(a±(t)a±(t + T)) ss = (a 2 ± {t)) ss e- x * T , 



<a±(t)a£(t + T))„ = 0, 



(<(*)<(< + r)) M = 2{M ± JV)5(r). 
In view of these results, the squeezing spectrum takes the form 

S£*(u) = 1 ± 2 K (a 2 ± {t)) ss Re / e^-^W =p 4k(M ± JV)i2e / 

Jo Jo 

/>oo 

±4(M±iV)i?e / e^ T 5(r)dT. 
Jo 

Thus upon performing the integration and using (|41[) , we easily find 



-(Azc — iujr) 



dr 



+ w 

Employing Eqs. (|30[) and (fTa[) -(|Tg|), the squeezing spectrum can be written as 

nout ( \ 2r 
£> + \UJ) — e 



1 + 



+ (l+/3S)(4+^)J "I" 2(l+/3 i )(l+^/4) 



and 



1 - 



Q„r r I - 4 ( 4 ./ 3 +/ 33 ).. l 3kA{3 2 e 2r 

AK \- £ + (l+^)(4+^)J 2(1+/3^)(1+/3V4) 



(l+/3*)(4+/3*). 



(61) 



(62) 



(63) 

(64) 

(65) 
(66) 
(67) 
(68) 



(69) 



(70) 



(71a) 



(71b) 
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so that on account of (|43|) the squeezing spectrum reduces at threshold to 



-,2r 



1 



K,A{A+f3 2 )e- 2r 
2(1+/3^)(1+/3V4) 



and 



3kA/3 



3K,Af3 2 e 2 



1 



(l+/3 2 )(l+/3 2 /4) 2(l+ / 3 2 )(l+ / 3 2 /4) 



(72a) 



(72b) 



L' 1 1 (l+^)(2+/3V2)^ 1 ~ 

In Fig. 8, we plot the squeezing spectrum [Eq. (|72b[) ] versus oj and (3. This plot indicates that perfect squeezing is 
attainable for oj — 0, f3 = and for any values of A, r, and k. To see the effect of the squeezed vacuum reservoir 
on the squeezing spectrum, we plot Eq. (|72b[) versus j3 for oj = 0, k = 0.8, A = 100, and for different values of the 
squeeze parameter r. It is easy to see from Fig. 9 that the presence of the squeezed vacuum reservoir increases the 
degree of squeezing significantly for [3^0. 

It is interesting to consider once more the case for which the nonlinear crystal is removed from the laser cavity. 
Thus upon setting e = A/i = (with [i ^ 0) in Eqs. (|71ap and (|71b|) , we obtain the squeezing spectrum to be 



and 



Sl ut {u) 



S^_ ut (oj) = e 



O r A(-2P+P*) i 



KA(4+l3 2 )e- 2r 
2(l+/3 2 )(l+/3 2 /4) 



2n[ Tl 



AjAfi+fi 3 ) 



3K.Af3 2 e 2 



(l +/ 32)( 4+(3 2 ) j 2(l+/32)(i +/ 32/ 4 ) 



ll + 



A(i/3+l33) 



0\ 2 



(73a) 



(73b) 



(l+/3 2 )(4+/?2) 

Expressions (|73a|) and (|73bp represent the squeezing spectrum of a coherently driven three-level laser coupled to a 
squeezed vacuum reservoir. In Fig. 10, we plot Eq. (|73b|) versus (3 for w = 0, k = 0.8, A — 100, and for different 
values of the squeeze parameter r. These plots show that the degree of squeezing of the output mode increases with 
the squeeze parameter r and almost perfect squeezing occurs for small values /3. 



V. PHOTON STATISTICS 



We now proceed to calculate the mean photon number of the cavity and output modes. The mean photon number 
of the cavity mode can be written as 

n= {a*(t)a(t}}. (74) 
Employing (|32a|) and its complex conjugate and assuming that the cavity mode is initially in a vacuum state, we get 

(a*(t)a(t)) = (F*(t)F(t)) (75) 
It can be readily established using Eqs. (|33aj) and (|33b|) along with (|26| and (|27| that 

(F*(t)F(t)) = 2A ~™ + £ (1 - e- 2A - 4 ) + 2 ^+^- £ (l _ e -^«). (76) 



(77) 



On account of Eqs. (|7a|) - (|7h)) . (1301) . CSt an( i (ESJ), the mean photon number turns out to be 

- _ [2e + n(e^ !)](! + /3 2 )(1 + /j 2 /4) + A(2 (3 + /3 2 /2 + /3 3 /2) _ _ 2A 

" 4(l + /3 2 )(l + /3 2 /4)( K -2e) + 2A(2/3-/3 3 ) 1 6 j 

[2e + n{er^ - !)](! + /3 2 )(1 + /3 2 /4) + A(2/3 - 3/3 2 /2 + /3 3 /2) _ 2A+ 

4(1 + /3 2 )(l+/3 2 /4)( K + 2e) + 2,4(4/3 + /3 3 ) 1 6 J ' 

Fig. 11 shows that the presence of both the parametric amplifier and the squeezed vacuum reservoir increase the 
mean photon number significantly for relatively small values of (3. On the other hand, it can be readily established 
that the mean photon number of the output mode is expressible as Q 

n ut = Kfi + (1 - k)N, (78) 

where n is given by Eq. J77|) . Fig. 12 shows that the mean photon number of the cavity mode is greater than that 
of the output mode. 
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VI. POWER SPECTRUM 



We finally calculate the power spectrum of the cavity and output modes. 



A. Power spectrum of the cavity mode 

The power spectrum of a single-mode light is expressible in terms of c-number variables associated with the normal 
ordering as 

/>oo 

S(u) = 2Re {a*{t)a{t + T)) ss e iuJT dr. (79) 



The solution of Eq. (p9|) can also be written as 

a(t + t)= E + (r)a(t) + E-(r)a*(t) + F(t + r), (80a) 

in which 

E ± {T)= l -{e- x -^±e- x n (80b) 

and 

F(t + r) =F+{t + r)+F_(t + T), (81a) 

with 

F±(t + r) = \ J e- x ^- T '\r(t + r 1 ) ± f(t + r'))dr'. (81b) 

Upon multiplying Eq. (|80ap by a* (t) and taking the expectation value of the resulting expression, we have 

(a*(t)a(t + t)) ss = E+(a*(t)a(t)) ss + E_(a* 2 (t)) ss . (82) 

It can be readily established using Eqs. (|32a[) . (|32b[) . (|33ap . and (|33b[) together with the properties of the noise force 
f(t) that 

/ 2\ / *2 N 2A~2V + e 2A + 2V- 

(a) ss = (a ) ss = — — . (83) 

In view of Eqs. J7HD, (|80b[) . and (JH3D , expression takes the form 

{a - (t)a{t + T »_ = M_21±£ e -- + 2C^Z£.^. ,84, 

Applying this result in (|79[) and performing the integration, the power spectrum of the cavity mode is found to be 

A-V + % A + V 

Employing Eqs. (|7ap - (|7h[) and (|30p . the power spectrum can be written in the form 



^ = l a^ + ^TT^ l - (85) 



K IJ2.T i\ , e i A(4-2/3+/3 2 +/3 J ) re f„— 2r -i\ i_ 

q( \ _ 4^ 1 ^ t 2 i ' (l+ffl)(8+2^) 4^ ^ 2 t (l+ffl)(8+2ffl) , , 

1.2 fc ^ (l+/3 2 )(4+/3 2 )J 12 + fc ^ (l+/3 2 )(4+/3 2 )J +UJ 

We see that the power spectrum is a sum of two Lorentzians centered at zero frequency, with a halfwidth of \ — 

£ + (1+^(4+^2) an d § + £ + (1+^(4+ ,32) • We realize that halfwidth of the two Lorentzians does not depend on the 
squeeze parameter. However, the halfwidth depends on e which represents the parametric amplifier. In Fig. 13 we 
plot the power spectrum of the cavity mode given by Eq. (|86p versus w for different values of e. These plots show 
that the halfwidth decreases with e. When the value of e increases from 0.2 to 0.3 the halfwidth decreases from 0.80 
to 0.75. 
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B. Power spectrum of the output mode 



The power spectrum of the output mode is expressible in terms of c-number variables associated with the normal 
ordering as 



S out (cu) = 2Re / (a* out (t)a out (t + r)) ss e^ T dT, 
Jo 



where 



w(*) = y/Ka(t) - a in (t). 



Now with the help of Eq. ([88)1 , we can write 



(°Cmt(t) a out(t + t)) ss = K(a*(t)a(t + t)) ss - ^/n(a* (t)a in (t + r)) s 
- \/K(a* in (t)a(t + t)) ss + (a* n {t)a in (t + r)) f 



It can be readily verified that 



y/K{a*(t)oiin(t + t)) ss = 0, 



(87) 



(89) 
(90) 



V^(a* n (t)a(t + t)) ss = ~[{M + N)e- X - T + (M - N) e - X + T ], 



(91) 



Thus on account of Eqs. 



(4(*)<»*.(Hr))„=%). 
0, (jniD, and ((§211, Eq. ([Ml) takes the form 
~2(A-V)+£ M + N 



(a* out (t)a ut(t + t)) ss = k 



4A_ 2 
2{A + V) - e M-N 



-X-T 



-X+T 



N5(t). 



4A+ 2 

On substituting this result into (|8T[) and carrying out the integration, we get 

_ A-V + §-X-(M + N)} K [A + V-I-X + (M-N)} 

With the aid of Eqs. (|30[) and (f?aj) - (|7h[) , we find the power spectrum of the output mode to be 



N. 



r,e _ A(2J3-t3 3 ) s 2 r , A(±+f3 2 ) 1 
h [^2 (l+/3 2 )(8+2 / 3 2 )V e T (l+/3 2 )(8+2/3 2 )] 

1.2 6 + (l+/3 2 )(4+/3 2 )J + W 



4-(f + 



A(4ff+j9 J ) y-2r 



2 ^ (l +/3 2)( 8+ 2/32), 



3A/3 2 ] 
(l+/3 2 )(8+2/3 2 )J 



1.2 ^ fc ^ (l+/3 2 )(4+/3 2 )J 



sinh 2 r. 



(92) 



(93) 



(94) 



(95) 



Expression (|95[) indicates that the power spectrum consists of two Lorentzians and a flat spectrum. Equations. (|86[) 
and (|95|) show that the width of the Lorentzians in the power spectrum of the cavity and the output modes are the 
same. 



VII. CONCLUSION 



We have considered a degenerate three-level laser containing a parametric amplifier and coupled to a squeezed vac- 
uum reservoir. Applying the pertinent master equation, we have obtained stochastic differential equations associated 
with the normal ordering. Using the solutions of these equations, we have calculated the quadrature variance, the 
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mean photon number, and the power spectrum of the cavity and output modes of the system under consideration. 
We have also determined the squeezing spectrum of the output mode. 

We have seen that the effect of the squeezed vacuum reservoir and the parametric amplifier is to increase the 
mean photon number and the degree of squeezing of the cavity and output modes. We have also found that the 
squeezed vacuum reservoir increases the degree of squeezing significantly over and above the squeezing attainable 
from a degenerate three-level laser with a parametric amplifier [2| ■ It turns out that the degree of squeezing of the 
cavity mode is greater than that of the output mode for certain values of (3. On the other hand, the plots of the 
squeezing spectrum at threshold show that there is perfect squeezing of the output mode for = uj = and for any 
values of A, k, and the squeeze parameter r. We have also found that the mean photon number of the cavity mode 
is greater that of the output mode. Moreover, we have seen that the presence of the parametric amplifier leads to a 
decrease in the width of the power spectrum while the squeezed vacuum reservoir has no effect on the width. 
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FIG. 1: A degenerate three-level laser with a parametric amplifier and a squeezed vacuum reservoir. 
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FIG. 2: A plot of the quadrature variance [Eq. (I44b|l ] versus f3 and r for k, — 0.8 and A = 100. 




FIG. 3: Plots of the quadrature variance versus P for k = 0.8 and A = 100, for r = and e = [Eq. (|45b|) ] (dashed curve), for 
r — in the presence of the parametric amplifier at threshold [Eq. (I44b|l ] (solid curve), and for r = 1.0 and in the presence of 
the parametric amplifier at threshold [Eq. (|44b[) ] (dotted curve). 




FIG. 4: A plot of the quadrature variance [Eq. (|45b|) ] versus r and /3 for k = 0.8 and A=100. 



14 



0.9 
0.8 
0.7 



0.6 




0.05 0.1 0.15 0.2 



FIG. 5: Plots of the quadrature variance [Eq. (|45b[l ] versus (3 for n — 0.8 and A — 100 in the absence of the squeezed vacuum 
reservoir (solid curve) and in the presence of the squeezed vacuum reservoir with r = 1.0 (dashed curve). 
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FIG. 6: (a) A plot of the quadrature variance of the output mode [Eq. ()59b[l ] versus f3 for k — 0.8, A — 100, and r = 1. (b) A 
plot of quadrature variance of the cavity mode [Eq. ()44b[l ] versus /3 for k = 0.8, A=100, and r = 1. 
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FIG. 7: (a) A plot of the quadrature variance of the output mode [Eq. (|60b[) ] versus /3 for k = 0.8, A — 100, and r = 1. (b) A 
plot of quadrature variance of the cavity mode [Eq. (I45b[l ] versus (3 for n — 0.8, A=100, and r = 1. 
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FIG. 10: Plots of the squeezing spectrum [ (|73b|l ] versus j3 for ui = 0, k — 0.8, A = 100, and for different values of the squeeze 
parameter r. 
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FIG. 11: Plots of the mean photon number [Q77f] at steady state versus (3 for k — 0.8, A — 25 and (a) for r = and e = A/i = 
with jj, ^ 0, (b) for r — and e = 0.3, (c) for r = 1.0 and e = 0.3. 




FIG. 12: (a) A plot of the mean photon number [l|78p] at steady state versus /3 for k = 0.8, A = 25, r = 1.0, and for e = 0.3. 
(b) A plot of the mean photon number [ (|77[) ] at steady state versus /3 for k = 0.8, A — 25, r = 1.0, and for e = 0.3. 




FIG. 13: Plots of the power spectrum [Eg . (|86f) ] versus ui for A=100, j3 = 0.01, r — 1.0, k — 0.8, and for different values of e. 



